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Abstract — Dynamic stability of an AUV determines how the 
vehicle behaves when disturbed while initially travelling on a 
straight course with no control input. The hydrodynamic 
analysis is conducted on an AUV to predict the effects of 
viscosity on dynamic stability in CFD studies comparing the 
degree of stability at various velocities. The hydrodynamic 
coefficients are calculated using numerical approach of 
Computational Fluid Dynamics (CFD) study using ANSYS 
Fluent 15.0 a non commercial software. The CFD studies which 
have to be conducted using viscosity effects require high 
performance computers and are slightly labor intensive (e.g. k- 
E and k-o) models) than inviscid model. By avoiding the usage 
of viscous models the CFD analysis can be completed quite 
efficiently with minimum computational requirements. The 
current paper investigates the need of whether to consider the 
effects of viscosity or not in order to produce the required 
analysis in a shorter time. So the main aim of the paper is to 
demonstrate two different model tests numerically, drift angle 
and rotating arm tests through which a stability margin is 
calculated to find out whether the AUV is directionally stable or 
not. The tests are repeated at various velocities in both viscous 
as well as no-viscous models, so as to clearly note down the 
difference. By the end of the paper, we would try and suggest as 
to in which operating velocities (whether high or low) it is 
necessary to consider the viscous effects. 


Index Terms —Dynamic stability, CFD, drift angle test, 
rotating arm test. 


I. INTRODUCTION 

Autonomous underwater vehicles have been used for a 
limited number of tasks dictated by technology available. 
With the development of more advanced processing 
capabilities and high yield power supplies, AUV’s are now 
being used for more and more tasks, with roles and missions 
constantly evolving. The oil and gas industry uses AUV’s to 
make detailed maps of the seafloor before they start building 
subsea infrastructure, pipelines and subsea completions can 
be in the most effective manner with minimum disruption to 
the environment. 



Fig 1. Autonomous underwater vehicle 
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Scientists use AUV’s to study lakes, the oceans, and the 
ocean floor. A variety of sensors can be affixed to AUV’s to 
measure the concentration of various elements or 
compounds, the absorption or reflection of light, and the 
presence of microscopic life. Military applications include 
intelligence, surveillance, reconnaissance, mine 
countermeasures, anti submarine warfare, inspection and 
identification (search & rescue), navigational network nodes 
and payload delivery in critical times of warzone. 

The dynamic stability of an AUV determines how the 
vehicle behaves when disturbed while initially travelling on a 
straight course with no control plane input. Using only linear 
terms, solutions to the sway and yaw equations provide linear 
transfer functions permitting the review of the stability of 
motion. The condition for stability is simply 


To know the hydrodynamic coefficients Y v , N v , Y r and N r 
there are various methods like experimental, computational 
methods. Although experimental methods are quite accurate, 
it becomes a second choice for the designers with the advent 
of computers because experimental methods are costly and 
time consuming whereas computational results are reliable 
and can be easily achieved. [1] 

In this work the dynamic stability of AUV is found out 
using two different viscous models (k-E and inviscid )at 
various speeds to find the appropriate speeds at which k-E 
should be used so as to reduce computational time of the 
analysis. 

II. OVERVIEW 
A. Model Description 

The model used here is a prolate spheroid. The curvature of 
profile is taken in 2-d form and it is rotated about the 
centerline axis in AutoCAD to obtain the 3-d form. Length 
(L) of model is 2.04m and the diameter(D) is 0.4 m. Fins have 
been attached at the aft to simulate realistic flow around the 
AUV. Fin geometry is obtained by extruded aerofoil section 
NACA (0019). The span of fin is 0.4 m. 



Fig .1. AUV model 
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B. Computational domain (Boundary Conditions) 

The important aspect while generation of domain is to 
ensure that the walls have little effect on disturbance of flow 
around the hull. 

For DRIFT ANGLE TESTS: A cuboidal domain whose 
dimensions are 60 m x 24 m x 24 m respectively was created 
in ANSYS design modeler. Front face, upper face, lower face 
and side faces are constrained at 12xL distance from model 
while back face is constrained at 25xL distance from model. 
The head of AUV hull is placed at 10 m from the inlet. [5] 



Fig .2. Meshed domain for drift angle tests 


For ROTATING ARM TESTS: Flow domain is defined as a 
semi-cylinder with a hole at the center. The radius of hole is 
equal to L, cylinder radius is equal to 8xL and height of 
cylinder is equal to 3xL.The AUV hull is placed at varying 
radius so as to alter the angular velocity. 


JL 

Fig .3. Meshed domain for rotating arm tests 


Table 1.boundary conditions of computational domain 


Boundary 

Type 

Inlet 

Velocity inlet 

Outlet 

Outflow 

AUV hull 

No slip wall 

Remaining sides 

No slip wall 



D. Solver model description 

All the analysis is done using Ansys FLUENT 15.0 
software. In this study, straight-line towing and Rotating Arm 
tests are simulated by using CFD methods. FLUENT code is 
used as flow solver. FLUENT solves mass and momentum 
conservation equations for all flow problems. In this study 
fluid is water so that it is assumed that flow is incompressible 
and there is no heat transfer. [2] 

All the tests are performed using two models, Standard k-E 
model and in viscid model. The incompressible, isothermal 
Reynolds Averaged Navier Stokes (RANS) equations are 
solved in order to determine the Cartesian flow field 
(Ui =u,v,w) and pressure (p) of the water around the AUV 
hull 


( 2 ) 

dx t oxj pdxj \dxj dXiJ J dxj JL 


By time averaging the Navier Stokes equations to generate 
the RANS equations, 6 further unknowns have been created, 
termed the Reynolds stresses: 


Various turbulence models are used to provide solutions to 
the Reynolds Stresses in terms of known quantities to allow 
closure of the RANS equations. The k- E model is a 
commonly used turbulence model for engineering 
simulations due to its robustness and application to a wide 
range of flows. Along with these, inviscid tests are conducted 
using inviscid flow model, which involves removal of 
viscous terms from the governing equations. 


III. TEST PROCEDURE 

A. Drift Angle Tests 

In actual practice the model is towed at various angles of 
attack, measuring sway forces and yaw moment. The sway 
velocity is = U sin (I , where U is the steady state inlet 
velocity and /? is the angle of attack. The sway forces and 
moments are plotted against the sway velocity, the slope of 
curve at zero angle determines Y v and N v respectively. 
Similar procedure is followed by varying the angles of attack 
of AUV model (i.e, 0,3,6,9 and 12 degrees) about the Z axis. 
For each inlet velocities, k-E and inviscid models are used. 
The tests are conducted at 0.5m/s, 2 m/s and 5 m/s. 


C. Meshing 

Ansys Mesh tool is used for the meshing of computational 
domain. Non Uniform unstructured mesh elements have been 
used to divide the domain into small control volumes. Mesh 
structure has been designed denser in the vicinity of the 
model and gradually the elements size increases as it moves 
away from it. Such an approach enables simulations run 
economically. The analysis is done using Ansys Fluent 15.0 
software. Flow domain is defined large enough to minimize 
flow effects between model and boundaries. 
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Fig .4.Flow visualization in drift test 
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B. Rotating Arm Tests 

The Rotating Arm Technique measures the rotary 
derivatives, or more specifically the Yr and Nr on the models, 
which is a specially arranged towing tank and apparatus. The 
model is being fixed to the end of a radial arm and an angular 
velocity is applied on it by rotating the arm through a vertical 
axis fixed in the tank. The orientation of the model is such 
that its x-axis and z axis is kept normal to the radial arm and is 
attached to the arm at the model’s mid-length. As a result of 
this orientation, when the model revolves about the axis of 
the tank, it rotates at the rate r. At this point however, the 
transverse velocity component v of the model is always zero 
(yaw angle of attack (3=0) and the axial velocity component 
ul is same to its linear speed. The model is constantly rotated 
at a constant linear speed and various radii R. [4] 


The roots of this quadratic equation should be negative for 
controls fixed stability; C is considered the discriminant of 
dynamic stability. The basic criterion for dynamic course 
stability can be expressed as 


N r 

M 



( 8 ) 


V. RESULTS AND DISCUSSIONS 
A. COMPARISON 

In this study, drift angle tests and rotating arm tests are 
performed at various inlet velocities 0.5,2 and 5 m/s using 
k-E and in viscid models. 


Experimental rotating arm tests have some limitations, i.e., 
they require large specialized facilities. In order to determine 
the values of Y r and N r as r^O the radius(R) should be 
large in relation to the AUV length (L) and the model must be 
accelerated and tests performed within a single revolution to 
ensure the vessel is not disturbed by its own wash, this limits 
the duration of each run. By performing virtual tests in a 
numerical towing tank these limitations can be overcome. 


Table 1. Stability derivatives using inviscid model 


Velocity (m/s) 

Yv 

N v 

Yr 

N r ' 

0.5 

-0.0315 

0.0505 

0.0559 

0.1001 

2 

-0.0331 

0.0515 

0.1429 

0.1003 

5 

-0.0314 

0.0671 

0.1609 

0.0980 



Fig .5. Flow visualization in rotating arm test 

IV. CALCULATING DYNAMIC STABILITY 
The numerical tests conducted are by fixing the controls, i.e., 
open loop stability which is studied only in horizontal plane. 
Only the linear terms are sufficient to review stability of 
motion. For the vessel to be stable, the horizontal velocity 
component v and yaw angular velocity component r should 
approach to zero as time increases to large value. [3] 


= W x e a ^ + V 2 e a2t 

(5) 

= R x e a ^ + R 2 e a2t 

(6) 


The simultaneous solution of these two equations for v and 
r' yields a second order differential equation which leads to 
dynamic stability criteria, e = 2.718, V 1 ,V 2 ,R 1 and R 2 are 
constants of integration; o x and o 2 are stability indices. The 
relationship between the stability indices and stability 
derivatives is 

Act 2 + Bo + C = 0 
Where A = n z Ay 

B = -n z Y' v -A y N' r (7) 

c = Y' v N' r - (Y' r -A')N' V 


Table 2. Stability derivatives using k-E model 


Velocity (m/s) 

Y v ' 

N v 

Yr 

Nr 

0.5 

- 0.0243 

0.0519 

0.1954 

0.0919 

2 

- 0.0225 

0.0503 

0.1737 

0.1164 

5 

- 0.0337 

0.0563 

0.1780 

0.1197 


Tables 1 and 2 clearly demonstrate the differences in the 
non-dimensional stability derivatives. Upon substituting 
these values in the stability criterion, we can predict whether 
the vessel has initial dynamic course stability or not 

Table 3. Dynamic stability at various velocities 


Velocity (m/s) 

k-E 

inviscid 

0.5 

stable 

unstable 

2 

stable 

stable 

5 

stable 

stable 


It is clear from above results that there is an effect of viscosity 
on dynamic stability at low speeds, whereas at high speeds of 
operation there is no significant effect. 


B. Validation 


The CFD analysis performed in this study is validated with 
the IITC formula (1957) for skin frictional resistance. 


_ 0.075 

F log(Rn-2) 2 


(9) 


The model is placed at an angle of attack of 0°,and inlet 
velocities are 0.5, 2 and 5 m/s. 
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VELOCITY (m/s) 

Fig .6. Drag force vs. Velocity 
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From the above graph it can be shown that, the methodology 
followed in this study is sufficiently accurate, because the 
difference is negligible. 


VI. Conclusion 

In this paper it is shown that there is an effect of viscosity on 
dynamic stability at low speeds, but the effect at high speeds 
is not considerable. So it is advisable to use viscous models 
while analyzing flow at low speeds. By using inviscid model 
at high speeds there is a considerable reduction in simulation 
time which is highly advantageous in low specification 
computers 
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